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ABSTRACT 

We explore fallback accretion onto newly born magnetars during the supernova of massive stars. 
Strong magnetic fields (~ 10 15 G) and short spin periods (~ 1 — 10 ms) have an important influence 
on how the magnetar interacts with the infalling material. At long spin periods, weak magnetic fields, 
and high accretion rates, sufficient material is accreted to form a black hole, as is commonly found for 
massive progenitor stars. When B < 5 x 10 14 G, accretion causes the magnetar to spin sufficiently 
rapidly to deform triaxially and produce gravitational waves, but only for rs 50 — 200 s until it collapses 
to a black hole. Conversely, at short spin periods, strong magnetic fields, and low accretion rates, 
the magnetar is in the "propeller regime" and avoids becoming a black hole by expelling incoming 
material. This process spins down the magnetar, so that gravitational waves are only expected if the 
initial protoneutron star is spinning rapidly. Even when the magnetar survives, it accretes at least 
w 0.3M Q , so we expect magnetars born within these types of environments to be more massive than 
the IAMq typically associated with neutron stars. The propeller mechanism converts the ~ 10 52 ergs 
of spin energy in the magnetar into the kinetic energy of an outflow, which shock heats the outgoing 
supernova ejecta during the first ~ 10 — 30 s. For a small ~ 5Mq hydrogen-poor envelope, this energy 
creates a brighter, faster evolving supernova with high ejecta velocities ~ (1 — 3) x 10 4 km s _1 and may 
appear as a broad-lined Type Ib/c supernova. For a large > IOMq hydrogen-rich envelope, the result 
is a bright Type IIP supernova with a plateau luminosity of > 10 43 ergs s" 1 lasting for a timescale of 
- 60 - 80 days. 

Subject headings: gravitational waves — stars: magnetic fields — stars: neutron — supernovae: general 



1. INTRODUCTION 

"Magnetars" are a subset of neutron stars with dipole 
' magnetic fields as strong as B ~ 10 14 — 10 15 G (Dun- 
; can & Thompson 1992; Thompson & Duncan 1993). Al- 
though at an age of 1, 000 — 10, 000 years they have spin 
■ periods of P = 5 — 12 s, as measured from soft gamma- 
| ray repeaters and anomalous X-ray pulsars, it is an out- 
f^S . standing question of how rapidly they rotate when first 
born. Short initial spin periods (Pq ~ 1 — 10 ms) have 
y—i ' been favored theoretically so that the dynamo process 
kJ . that create s these strong magnetic fields may operate 
• £h ' efficiently dDuncan fe Thompson! Il992t lAkivama et al.l 
K> ; 120031: iThompson et all 120051 ). Motivated by this, many 
groups have investigated the possible impact of the spin- 
el ■ down of this newly formed magnetar in powering an ex- 
plosion (see, for example, Bodenheimer & Ostriker 1974; 
Wheeler et al. 2000; Thompson et al. 2004; Burrows 
et al. 2007; Dessart et al. 2008). Such short spin pe- 
riods may also be a source of ultra-high energy cosmic- 
rays (Arons 2003), create a colli mated relativistic flow as 
needed for gamma-ray bu rsts (jUzdensky fc MacFadvenl 
l2007t iMetzger et all 120101 and references therein), or 
produce a luminous supernova (Kasen & Bildsten 2010; 
Woosley 2010). 

An assumption of all of these studies is that the su- 
pernova which gave birth to the magnetar was successful 
in ejecting the majority of the progenitor star's enve- 
lope. This is clearly correct in many cases, since we 
know that neutron stars with more modest magnetic 
fields (~ 10 12 G) are created in supernovae. But it is 
possible that some subset of supernovae which produce 



neutron stars have small injected explosion energies. As 
is expected for massive stars that give rise to black holes, 
these would not be successful in ejecting the majority of 
the envelope and a sizable amount of fallback would oc - 
cur (as found for ss 25— 40M Q stars bv lHeger et all2003f ). 
In addition, even in cases where the majority of the en- 
velope is ejected, asymmetries in the explosion may still 
result in significant fallback. For these reasons, it is plau- 
sible that there exists a population of massive stars that 
give birth to magnetars that are subsequently subject to 
accretion of the envelope material. 

Another motivation for studying fallback accretion 
onto magnetars is the presence of magnetars near clus- 
ters of massive stars. SGR 1806 - 20 and CXOU 
J164710.2 — 455216 are associated with the clusters CI 
1806 — 20 and Westerlund 1, respectively, and are in- 
ferred to have had progenitor masses of « 40M Q (Figer 
et al. 2005; Bibby et al. 2008; Muno et al. 2006). Fur- 
thermore, the expanding H I shell around the magnetar 
IE 1048.1 - 5937 also argues for a « 30 - 4OM pro- 
genitor (Gaensler et al. 2005). Such massi ve stars are 
typically assumed to give rise to black holes (|Frverlll999l : 
iHeger et all 120031 ). although we note that this will de- 
pend sensitively on the details of mass loss during stel- 
lar evolution (Smith et al 2010; O'Connor & Ott 2011) 
and on whether these magnetars have binary progenitors 
(Belczynski & Taam 2008). It is therefore worth explor- 
ing whether the presence of a highly-magnetized neutron 
star qualitatively changes the outcome of the collapse of 
massive stars. 

In the following study we explore the interaction of 
newly born magnetars with supernova fallback. We be- 
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gin in ^ by discussing the parameter space in which 
we expect fallback to be important. In Sj3]we calculate 
the time-dependent spin evolution of these magnetars. 
These results are used in $4] to explore whether a newly 
formed magnetar accretes sufficient material to become a 
black hole, as a function of the initial spin, magnetic field, 
and amplitude of the fallback accretion. We also discuss 
whether these magnetars will be spinning rapidly enough 
to produce gravitational waves via triaxial instabilities. 
In Sj5] we show that material expelled in the propeller 
regime collides with outgoing supernova ejecta, creating 
a more powerful supernova. We conclude in Ej6] with a 
summary of our results. In the Appendix we explore the 
physics of neutrino-cooled accretion columns onto mag- 
netars. 

2. FALLBACK VERSUS OUTFLOW 

Before we investigate the effects of fallback accretion, 
it is pertinent to discuss when fallback is expected. Al- 
though these arguments are strictly applicable for only 
one-dimension, and we expect a multi-dimensional flow 
to provide more opportunities for fallback, this gives 
some intuition about how fallback depends on the ac- 
cretion rate, spin, and magnetic field strength. 

As the rapidly rotating, newly born magnetar spins 
down, it goes through stages in which it emits en- 
ergy in dipole spindown radiation and a neutrino-driven, 
magnetically dominated wind (Thompson, Chang, & 
Quataert 2004), both of which may hinder accretion. For 
a magnetar with a dipole magnetic moment /x and spin 
fi, the spindown luminosity is 



Ldip — 



6c 3 

= /i/10 33 



= 9.6 x 10 4 V^-Pf 4 er g s s ~\ 



(1) 



where /i33 = /x/10 G cm , as is appropriate for 
a neutron star with a 10 15 G magnetic field, and 
P = 2tt/Q, = 1P\ ms. Assuming this luminosity is car- 
ried by a relativistic wind, the associated pressure at a 
radius r is pdip = £dip/47rcr 2 . Fallback accretion exerts 
an inward ram pressure, and for the case of spherically 
symmetric accretion at a rate M onto a mass M, this is 
given by 

M (2GM\ lj2 

Pram = — ( — — ) ■ (2) 



Since p dip 



and p r 



x 



r 5 / 2 , the spindown lumi- 



nosity always wins at sufficiently large radii. If the fall- 
back accretion is already proceeding and then the spin- 
down luminosity is to disrupt this accretion flow, we can 
ask what is the critical accretion rate above which the 
fallback ram pressure dominates at the magnetar radius 
R. This gives 
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where M u = M/lAM e and R 12 = R/12 km. This 
accretion rate is well-exceeded in all cases we consider. 

During the Kelvin-Helmholtz cooling epoch for the 
newly born magnetar, deleptonization and thermal neu- 
trino losses create a neutrino-driven wind that is magnet- 



0.01 



0.001 = 



• s 0.0001 i 




P (ms) 

Fig. 1. — The critical accretion rate, above which fallback dom- 
inates, as a function of the spin period. We consider two physical 
processes for inhibiting the fallback: dipole spindown radiation 
(denoted by Mdj p>cr it and given by eq. [3]), and a neutrino-driven 
wind (denoted by Af uc rit an d given by eq. [5]). In each case we 
vary the radius by a factor of 100 (as shown by the shaded regions) 
to represent uncertainty in the radius at which the accretion flow 
first comes into contact with this outgoing energy. 

ically flung by the magnetar's dipole field. For a mass 
loss rate M u , the luminosity that goes into this process 
is (Thompson et al. 2004) 



L„ = 



2/5 



= 4.5 x 10 5 V /5 



p-8/5 ,V3/5 
r \ IVI v,-3 



ergs s 



(4) 



where M v ,-3 = Af„/10~ 3 M Q s -1 . Repeating the above 
analysis of assuming this is a relativistic wind and com- 
paring to the ram pressure at the magnetar surface, we 
derive a critical accretion rate 



• 2M„ ( R \ 



1/2 



2/5 



= 8.6 x 10-Vaf Pi ^M^M^'^R^ M s~\ 

(5) 

This limit is a little more stringent than the one derived 
for dipole spindown (eq. [3]). Indeed some of the lower 
fallback rates we consider are exceeded by this. When 
Thompson et al. (2004) follow the spindown from a 
neutrino-driven wind, they find modest amounts of spin- 
down (an increase in the spin period of ~ 5 ms) even 
for the most extreme conditions. If there is a phase of 
spindown from this, it just amounts to different initial 
conditions from the perspective of our study. Thus, we 
neglect these effects in our time-dependent spin calcula- 
tions. 

In Figure[T]we summarize the parameter space in which 
we expect fallback to be important. This shows that the 
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fallback ram pressure dominates for accretion rates above 
~ — lO~ 2 Af s -1 , depending on the process that 

is inhibiting the fallback. Comparing with the fallback 
found in numerical studies by MacFadyen et al. (2001) 
or Zhang et al. (2008), this implies a massive progeni- 
tor (in the range of ~ 20 — 4OM for solar metallicity 
and the progenitor models of Woosley et al. 2002) and 
a low explosion energy (higher explosion energies lead 
to weaker fallback; Dessart, Livne, & Waldman 2010). 
Although it is not well-known how progenitor mass and 
explosion energy correlate with magnetar creation, even 
with these limitations, there is a wide parameter space 
where fallback onto a magnetar seems inevitable. 

Even in cases where these scalings appear to argue that 
fallback is inhibited, it is still worthwhile to investigate 
fallback on account that (1) the neutrino-driven wind 
only lasts ~ 10 s while the fallback occurs on a > 1000 s 
timescale (reflecting the dynamical time of the progeni- 
tor) and (2) the neutrino-driven wind is highly asymmet- 
ric. Therefore, even if the wind excavates some region of 
the progenitor, there is ample opportunity for fallback at 
other angles. We thus expect that in higher dimensions 
the strength of fallback is typically greater than what we 
assume for our one-dimensional arguments. 

3. SPIN EVOLUTION DUE TO FALLBACK ACCRETION 
3.1. Accretion Versus Expulsion 

The initial spin period of newly-born neutron stars de- 
pends on both the spin profile of the progenitor star and 
subsequent processes that add, subtract, and redistribute 
angular momentum. Fryer & Heger (2000) performed 
smoothed particle hydrodynamics simulations using a ro- 
tating progenitor model from Heger et al. (2000) and 
estimated an initial protoneutron star (PNS) spin period 
on the order of 100 ms. It is, however, not clear how they 
defined the extent of the PNS (see discussion in Ott et 
al. 2006). The subsequent cooling and contraction to a 
radius of ~ 12 km resulted in Po ~ 2 ms. Fryer & Warren 
(2004) subsequently estimated neutron star spin periods 
by assuming that the angular momentum of the inner 
IMq is conserved as the PNS cools and contracts to a 
neutron star, finding periods of ~ 1 — 17 ms depending on 
the progenitor model. Thompson et al. (2005) studied 
the action of viscous processes in dissipating the strong 
rotational shear profile produced by core collapse in a 
range of progenitors and for different initial iron core pe- 
riods. They showed that for rapidly rotating cores with 
postbounce periods of < 4 ms, viscosity (presumably due 
to magnetic torques via the magnetorotational instability 
or magnetoconvection) spins down the rapidly rotating 
PNSs by a factor of ~ 2 — 3 in the early postbounce 
epoch. Ott et al. (2006) systematically studied the con- 
nection between progenitors and final neutron star spin, 
generally finding Po ~ 0-5 — 10 ms and solid body ro- 
tation in the PNS core for progenitors with precollapse 
periods < 50 s. We therefore consider initial magnetar 
spin periods in this range for our present study. 

Subsequent to the initial spin period being set as de- 
scribed above, the neutron star may be subject to fall- 
back accretion. Accretion comes under the strong in- 
fluence of the star's dipole field at the nominal Alfven 
radius r m — ^ 4 / 7 (GM) _1/,7 M -2 / 7 , where fi is the dipole 
magnetic moment of the magnetar. For typical magnetar 



parameters 

r m = Ufi% 7 M~l /7 MZ* /7 km, (6) 

where M_2 = M /10 -2 M Q s _1 , and the prefactor to 
r m can vary depending on the details of the interac- 
tion between the flow and magnetic field (jGhosh fc Lambl 
119791: lAronsI 119861 11993ft . The other critical radius, 
set by the magnetar's spin Q, is the corotation radius 
r c = (GM/n 2 ) 1 / 3 , 

r c = UM^P^ 3 km. (7) 

Roughly speaking, one expects that for r rn < r c , material 
is funneled by the magnetar's dipole field before accreting 
onto the magnetar's surface, while when r m > r c , ma- 
terial must spin at a super-Keplerian rate to come into 
corotation with the magnetar and is thus expelled (the 
"propeller regime," Illarionov & Sunyaev 1975). Setting 
r m > Tc gives a critical accretion rate 

M < 6.0 x 10~ 3 f43 M il /3p i 7/3 M Q s_1 - ( 8 ) 

Comparing to the 25 M© collapsar models of MacFadyen 
et al. (2001), they find early-time accretion rates of 
10 -4 — 10~ 2 Mq s _1 by just varying the injected ex- 
plosion energy by (0.255 — 1.2) x 10 51 ergs. Whether a 
magnetar is in the propeller regime or not is therefore 
very sensitive to how energetic the supernova is. 

This simplistic picture is not the complete story, as 
has been detailed by a great many theoretical studies of 
accretion onto magnetic stars (see for example, Pringle & 
Rees 1972; Lynden-Bell & Pringle 1974; Ghosh & Lamb 
1979; Aly 1980; Wang 1987; Shu et al. 1994; Lovelace et 
al. 1995, 1999; Ikhsanov 2002; Rappaport et al. 2004; 
Ek§i et al. 2005; Kluzniak & Rappaport 2007; D'Angelo 
& Spruit 2010). More recently, numerical simulations 
have also been used to investigate this problem (Hayashi 
et al. 1996; Goodson et al. 1997; Miller & Stone 1997; 
Fendt & Elstner 2000; Matt et al. 2002; Romanova et al. 
2003, 2004, 2009). For our present work, we implement 
a simple model largely based on that used by Ek§i et 
al. (2005), as described below. Their prescription has 
the advantage of being applicable and continuous over 
a wide range of parameters, while capturing the main 
expected features of the propeller regime. 

In cases where r c > r m > R, the inflowing material 
is channeled onto the magnetar poles where it shocks 
and neutrino cools. We save a more detailed treatment 
of the physics of this process for the Appendix, since 
it does not have a direct bearing on our results for the 
time-dependent spin, which we consider next. 

3.2. Time- Dependent Spin From Fallback Accretion 

Given this picture of accretion and expulsion described 
above, we solve for spin evolution under the influence of 
fallback accretion by integrating the differential equation 

I^jT = ^dip + A^acc, (9) 

where / = 0.35Afi? 2 is the moment of inertia (Lattimcr 
& Prakash 2001), and N^ip and 7V acc are the torques from 
dipole emission and accretion, respectively. As discussed 
in £j2l we ignore spindown from neutrino-driven winds in 
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equation Q. The dipole spindown torque is given by 



Ndip = 



6c 3 



= -1.5 x 10 4 ^ 3 Pf J ergs. (10) 



We assume that the magnetar is rotating as a solid body, 
as is likely the case within ~ 1 s of colla pse since the MRI 
([Thompson et al.ll2005t lOtt et al 120061) or low -T/\W\ in- 
stabilities (jWatts et alJl2005t lOtt etaLl 12005) will limit 
differential rotation. When r m > R, material leaves the 
disk with the specific angular momentum at a radius r m . 
Depending on the relative positions of the Alfven and 
corotation radii, this can either spin up or spin down the 
magnetar, so we write the torque as 

N acc = n{Lu)(GMr m ) 1 ' 2 M if r m > R, (11) 

where n{uS) is the dimensionless torque which depends 
on the fastness parameter us = Q/ (GM/r^) 1 / 2 = 
(r m /fc) 3 ^ 2 ■ Ek§i et al. (2005) discuss different ways 
in which n{ui) can be set, but for simplicity we take 
n = 1 — u). This has the advantage that the torque goes 
to zero at the corotation radius, is continuous for all u>, 
and goes negative when r rn > r c , corresponding to the 
spin down which occurs during the propeller regime. As 
ui gets larger, this prescription gives increasingly strong 
spindown, consistent with the more detailed simulations 
of Romanova et al. (2004). When r m < R we set the 
torque to 

A^acc = (l-0/0 K )(GMi?) 1/2 M if r m <R, (12) 

where fix — (GM/R 3 ) 1 / 2 . The prefactor is included 
to ensure that torque is continuous for all values of r m . 
The disadvantage is that since the prefactor is < 1, it 
will underpredict the amount of torque, but this does 
not change our main conclusions, as we discuss in §4.2. 

As we integrate the spin in time, we keep track of 
the magnetar's rotation parameter, f3 = T/\W\, where 
T = in 2 /2. We use the prescription given in Lattimer 
& Prakash (2001) for \W\, 



\W\ » 0.6Mc 2 



GM/Rc 2 



l-0.5(GM/Rc 2 )' 



(13) 



We keep R fixed even as M changes, which is roughly 
consistent with most equations of state, except when M 
gets near its maximum value (Lattimer & Prakash 2001). 
When j3 — 0.5, the neutron star is at breakup and can- 
not accept further angular momentum. Even prior to 
this, dynamical bar-mode instabilities occur for j3 > 0.27 
(Chandrasekhar 1969), and secular instabilities for j3 > 
0.14, driven by gravitational radiation reaction or viscos- 
ity (Lai & Shapiro 1995). Since the dynamical bar- mode 
instability is guaranteed to radiate and/or hydrodynami- 
cally re-adjust angular momentum, we set iV acc = when 
(i > 0.27. We ignore changes in spin due to the secular 
instabilities since growth timescales are uncertain and 
may be suppressed by competiti on between viscosity or 
gravitational radiation reaction (|Lai fc Shapirol lT995) . 

We parameterize the fallback accretion rate to mimic 
the results of MacFadyen et al. (2001) and Zhang et al. 
(2008). This can roughly be broken into two parts. At 
early times it scales as 



Afcarly = TjltT^M, 



where rj » 0.1 — 10 is a factor that accounts for different 
explosion energies (a smaller r\ corresponds to a larger 
explosion energy), and t is measured in seconds. The late 
time accretion is roughly independent of the explosion 
energy and is set to be 



AL 



late 



= 50t" 5/3 M m s _1 . 



(15) 



The accretion rate at any given time is found from com- 
bining these two expressions 



M 



Jw carly 



/if -1 



(16) 



The mass of the neutron star increases at a rate M when 
I'm < f c and is set fixed when r m > r c . For comparison, 
we also integrate M for all values of r m to follow how 
much matter the magnetar would have accreted if not 
for the propeller mechanism. 

Equation (fT6|) reflects fallback of the envelope, but 
most likely this material must pass through a disk be- 
fore finally accreting onto the magnetar. To test this 
hypothesis and explore whether this leads to a quanti- 
tative change of the accretion rate, we built one-zone, 
a-disk models (similar to Metzger et al. 2008) using 
the angular momentum profiles of the massive, rotating 
progenitors of Woosley & Heger (2006) simulated with 
GRID (O'Connor & Ott 2010). Our general finding was 
that (1) there is sufficient angular momentum to form a 
disk, and (2) the disk is nearly steady-state, where the 
accretion rate onto the star differs from the infall rate by 
no more than a factor of ~ 5 (and this scales with the 
a-viscosity, with a larger a resulting in higher accretion 
rates), and (3) the radius of the disk is typically well 
outside of the Alfven radius. We therefore consider the 
mediation of the disk to be degenerate with rj and use 
the direct infall rates as described above. 

In Figure [2j we compare integrations of equation ([9]) 
for values of r\ = 0.1, 1, and 10. The top panel shows the 
accretion rate given by equation ([T6|) . The middle panel 
plots the time-dependent spin period. The bottom panel 
plots the fastness parameter, which reflects whether or 
not the magnetar is in the propeller regime. For r\ = 0.1, 
only 0.25 M is accreted out of a potential amount of 
accretion of 1.55 M Q , and for r\ = 1 only 1.03 M Q is ac- 
creted out of a potential amount of 3.15 M . Therefore 
both these cases are able to avoid becoming a black hole 
via the propeller mechanism (assuming a maximum neu- 
tron star mass of 2.5M©). In contrast, the 77 = 10 case 
(which corresponds to a lower-energy explosion) accretes 
3.45 M Q out of 6.41 M Q , which means it likely becomes 
a black hole. Since the accretion rate is highest at early 
times, black hole formation happens rather quickly dur- 
ing the runs, at « 34 s and « 46 s for maximum neutron 
star masses of 2.5A/ and 3M , respectively. 

In each of these cases, the spin eventually reaches an 
equilibrium value that simply tracks M with Set- 
ting r m = r c , we calculate an equilibrium spin period, 



P cq = 2V /7 (GAf)- 5 / 7 M" 3 / 7 



= 5.8/x 



6/7 -5/7 

33 1V1 IA 



M 



-3/7 



(17) 



s 



(14) 



where Af_ 4 = A'//10" 4 M q s" 1 . 

4. MAGNETAR VERSUS BLACK HOLE FORMATION 
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Fig. 2. — The spin evolution of a magnetar with B = 10 15 G and 
an initial spin period of Pq = 1 ms. We compare values of r\ = 
0.1, 1, and 10, demonstrating the strong effect early-time accretion 
can have. The top panel shows the time-dependent accretion rate, 
the center panel shows the spin period, and the bottom panel shows 
the fastness parameter u>, where uj > 1 corresponds to the propeller 
regime and w < 1 corresponds to accretion. 
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Fig. 3. — Contours show the amount of mass accreted (in so- 
lar masses) for different initial spin periods and magnetic field 
strengths (all for r\ = 1). For every case we assume an initial 
magnetar mass of 1.4Mq with a radius of 12 km. If the propeller 
regime did not expel material, then 3.15A/q would have been ac- 
creted. 



4.1. The Amount of Mass Accreted 

The example models in the previous section demon- 
strate that the amount of mass accreted by the magne- 
tar depends strongly on whether the propeller regime is 
reached. Therefore, whether or not a magnetar even- 
tually becomes a black hole depends on its initial spin 
period and magnetic field. This is in stark contrast to 
neutron stars with dynamically unimportant magnetic 
fields whose fates simply depend on the properties of 
the supernova and the compactness of th e stellar core 
(IZhang et all l200i IQ'Connor fe"Ottl I20TI . To explore 
these correlations, we plot contours for the amount of 
mass accreted as a function of the initial spin period and 
magnetic field in Figures [3] and [4] for values of 77 = 1 and 
0.1, respectively. In the r\ — 1 case, a magnetar remains 
for only a small fraction of the initial conditions (this of 
course depends on the value of the maximum neutron 
star mass). For 77 = 0.1, a magnetar is expected for the 
majority of the parameter space. Since these two values 
of 77 correspond t o a factor of ~ 2 differen ce in the initial 
explosion energy (jMacFadven et afll2001[ ) , these compar- 
isons demonstrate just how sensitive the outcome is to 
this quantity. 

The general trend is that at high magnetic fields and 
small periods, there is less mass accretion due to the 
propeller mechanism being stronger. Nevertheless, there 
are also some subtle differences from this trend that are 
due to the interaction of a time-dependent accretion rate 
with the changing spin. For example, in Figure [3] we 
see that the minimum accreted mass occurs near B » 
10 15 G and Pq « 0.7 ms, and the accreted mass actually 
increases for stronger magnetic fields, contrary to our 
intuition for when the propeller mechanism should be 
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Fig. 4. — The same as Fig. [3] but with 77 = 0.1. For this case, 
if the propeller regime did not expel material, then 1.55M0 would 
have been accreted. 

strongest. To explore what is happening here, we plot 
the spin evolution for a collection of different magnetic 
fields and initial spin parameters in Figure [SJ We can see 
that at sufficiently strong magnetic fields, the propeller 
is so strong that the star quickly spins down during the 
first w 20 s (dotted line). At this point the accretion rate 
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Fig. 5. — Time-evolution of the spin period and fastness param- 
eter ui for a diverse selection of models. The accretion rate cor- 
responds to r\ = 1 for all cases (the solid line from the top panel 
in Fig. [2{. The low magnetic field (dot-dashed lines) and slowly 
spinning (dashed line) cases exceed a mass of 2.5Mq at fti 180 s, 
at which point they most likely become black holes. See the text 
for further discussion of the features exhibited here. 

has increased dramatically, and the star now accretes and 
spins up until about w 200 s. It is due to this stage that 
the magnetar accretes more than was expected. 

One takeaway message of this parameter survey is that 
for all these models at least « 0.3Mq is accreted. There- 
fore magnetars that are subject to the conditions of be- 
ing born within a massive star should on average be more 
massive than the 1.4M Q typically associated with neu- 
tron stars. Measuring the masses of magnetars would 
therefore be useful for constraining whether some are in- 
deed born in massive progenitors. 

4.2. Prospects for Gravitational Wave Production 

These results also have bearing on whether a young 
magnetar should be expected to be an important gravi- 
tational wave source (as discussed in Corsi & Meszaros 
2009, and references therein). For this to occur, it must 
be spinning sufficiently quickly that dynamical bar-mode 
instabilities or secular instabilities are excited. To ex- 
plore this, we plot the spin parameter /3 for a selection of 
models in Figure |6] The majority of the parameter space 
we probe experiences some time in the propeller regime, 
spinning down the magnetar, and making gravitational 
wave emission unlikely. For magnetic fields < 5 x 10 14 G, 
the magnetar is spun up by accretion sufficiently to pro- 
duce gravitational waves, but the accretion then quickly 
leads to collapse to a black hole. This is seen in the 
bottom panel of Figure El where (3 > 0.14 for a time, 
but then exceeds a mass of 2.5M Q at « 180 s. This 
model never exceeds /3 = 0.27, but this is an artificial 
effect of the 1 — Q/CIk factor for the torque prescrip- 
tion (see eq. [10]). If we instead assume the magnetar 
accreted with the specific angular momentum at its sur- 
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Fig. 6. — The spin parameter for a selection of models, all with 
77 = 1. The dotted lines denote p = 0.14 (the critical value for 
secular instabilities) and j3 = 0.26 (the critical value for dynamical 
instabilities). Gravitational wave emission is only expected above 
the dotted lines, and is generally seen for extremely short initial 
spin periods (Po < 0.7 ms) or small magnetic field strength(B < 
5 x 10 14 G). 

face of (GMR) 1 / 2 , /3 = 0.27 would be easily reached. We 
estimate the timescale for gravitational wave emission by 
integrating the early time accretion law, 



i gw = 140)7 



-2/3 



Mn 



1.1M 



2/3 



(18) 



where M max is the maximum neutron star mass before 
black hole formation and Mq is the initial neutron star 
mass. The accretion peaks on a timescale 



(19) 



t p = 150rT 6/13 s, 

which is found by equating equations (Till and (fl"5)) . So 
our two conditions for equation (fTg)) to be valid are that 
B < 5 x 10 14 G and t gw < t p . If B > 5 x 10 14 G we 
don't expect appreciable spinup and gravitational wave 
emission, and if i gw > t p then the gravitational wave 
emission timescale is merely « t p . 

5. PROPELLER-POWERED SUPERNOVAE 

In cases that do not collapse into black holes, the 
material expelled by the propeller mechanism collides 
with the supernova ejecta. This shock-heats the enve- 
lope and increases the energy budget of the supernova. 
We next estimate the observable signature of such pow- 
ering. The process we describe here is decidedly different 
from what was explored by Kasen & Bildsten (2010) and 
Woosley (2010), who used dipole spindown luminosity 
to heat and power a more luminous supernova. In their 
case the dipole spindown takes place on sufficiently long 
timescales that it can directly power an extremely lumi- 
nous supernova. As we discuss below, the majority of the 
energy from the propeller-mechanism is injected during 
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the first ~ 10 — 30 s, so it can be treated as a sudden 
impulse of energy at early times. The majority of this 
energy is therefore lost to adiabatic expansion and not 
seen directly in the peak luminosity. Nevertheless, the 
energy can accelerate the supernova ejecta to high veloc- 
ities of up to ~ (1 — 3) x 10 4 km s -1 , which are observable 
in the spectra and alter the lightcurve shape. 

5.1. Propeller Energy Budget 

The expelled material carries a kinetic energy equal to 
the spindown energy of the magnetar. To help power the 
supernova, this material must climb out of the magne- 
tar's gravitational well, so we estimate the kinetic lumi- 
nosity of the propeller material as 



iprop = ~N acc Q - GMM/r m , 



(20) 

where the negative sign in the first term is because we 
have defined A acc to be negative when the magnetar is 
spinning down (eq. [11J ). With this equation we have 
assumed that the majority of the outflow originates from 
the inner edge of the disk. While this is a reasonable 
assumption, it also means that the material has to travel 
the furthest out of the potential well. If material can 
leave the disk at larger radii, it will require less energy 
to do so, thus this represents a lower limit. The total 
energy that can possibly be put into expelled material is 
limited by the magnetar rotation, 



E IC 



1 



-m 2 = 2.8 x 10 



12 



(21) 



In some cases the early-time accretion may even spin the 
magnetar up to sub-millisecond spin periods before the 
propeller mechanism begins. In these cases the magnetar 
stores the accretion energy in its spin, which is tapped via 
the propeller mechanism to help power the supernova. 

In Figures [Jj and [8] we quantify the luminosity of the 
propeller mechanism as well as what fraction of E mt is 
able to be tapped by this process. The top panels of 
each figure show L prop (eq. [20]), L^ip (eq. pQ), and the 
radioactive decay of 0.3M Q of 56 Ni as a function of time. 
The propeller-powering only lasts ~ 10 — 30 s until r rn ~ 
r c . At this point the magnetar is not spinning sufficiently 
rapidly to expel material to infinity and the luminosity 
quickly shuts off. The bottom panels of Figures [7J and [5] 
show the integrated energy as a function of time, 



Ei{t) = / Li(t)dt, 



(22) 



where i stands for either the propeller luminosity or 
dipole luminosity. In Figure [7J less than ~ 20% of the 
rotational energy goes into expelling material. For a 
stronger magnetic field the propeller regime is more ex- 
treme, and nearly all of the rotational energy is converted 
into energy of outflowing material, as shown in Figure 
In either case, this additional energy may be greater than 
the typical supernova energy E sn of ~ 10 51 ergs. In the 
following sections we explore how this additional energy 
alters the properties of the supernova depending on the 
mass of the envelope material, representative of Type 
Ib/c and Type IIP supernovae. 

5.2. Low-Mass Envelopes 
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Fig. 7. — The luminosity emitted by the magnetar, either from 
kinetic energy of magnetically expelled material L prop (solid lines), 
dipole radiation Ldip (dashed lines), and, for comparison, decay of 
0.3Mq of 56 Ni (dotted lines). The initial parameters are Po = 1 ms 
with r\ = 1 with a magnetic field of 10 15 G. The bottom panel 
shows the integrated energy as a function of time for each case. The 
total energy from radioactive heating is not sufficient to appear on 
the bottom panel. 




1000 10 4 
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Fig. 8. — The same as Fig. [7] but with a magnetic field of 
5 X 10 15 G. This example is more strongly in the propeller-regime 
at early times, so that a larger fraction of the spin energy goes into 
the kinetic energy of the outflow. 
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We consider a supernova with initial energy E sn , ejecta 
mass M e j, subject to a sudden impulse of energy E plop . 
As we will show, the observational impact of E pTop de- 
pends strongly on the ejecta mass and composition. In 
this section we focus on the properties of an event with a 
hydrogen-deficient envelope and a mass M e j < 5M@, as 
is expected for the progenitors of Type Ib/c supernovae 
that have lost a large fraction of their envelope (includ- 
ing all of their hydrogen) t o stellar winds, b inary mass 
transfer, and/or outbursts (jSmith et al.H20lTI ). 

The collision of the propeller material with the super- 
nova ejecta shock heats and accelerates the ejecta. For a 
total energy i? to t = E sn + E prop , the final velocity, with 
which it coasts for the remainder of the expansion, is 



v f 



(2E tot /M ej ) 1 / 2 = 2500£ 5 1 2 /2 5 Af- 1/2 km s 



(23) 



where E 52 . 5 = E tot /3 x 10 52 ergs and M 5 = M oj /5M Q . 
The diffusion timescale of photons from this hot, expand- 
ing material is given by 



13.78w/c 



1/2 



= U$*M! /A E£(* days, (24) 



where k is the opacity, which we scale to kq.i = 
k/0.1 cm 2 g^ 1 (the typical opacity used for a gray calcu- 
lation, Pinto & Eastman 2000), and the factor of 13.78 
comes from detailed analytic studies of Ty pe I super- 
novae (|ArnetH 119821 IPinto fc Eastman! [2000T) . The shell 
becomes optically thin on a timescale 



tr 



3 M, 



1/2 



47T V 



22^1 M 5 E^i 2 



days. (25) 



The diffusion approximation we will use is not applica- 
ble after this time. The increased velocity creates a faster 
and more luminous supernova, but this higher luminosity 
is not directly from energy input from the propeller mech- 
anism. Instead, since the explosion velocity is higher, the 
diffusion time (eq. [24]) is shorter, and the 56 Ni decay is 
being probed at earlier times. 

To understand the corresponding lightcurve created 
by propeller energy being injected into the explosion, 
we construct a simple, one-zone model of the expan- 
sion, cooling and emission, following the mathematical 
framework of Li & Paczyhski (1998; also see Kulkarni 
2005, Kasen & Bildsten 2010). For an expanding shell, 
the internal energy E- m t satisfies the differential equation 
(rewritten from eq. [9] of Li & Paczyhski 1998) 

\jt ^ Eint ^ = V°p(*) + £nuce- i/tT - L(t), (26) 



where 



L(t) = E- mt (t)t/4 



(27) 



is the emitted luminosity. The nuclear luminosity 
includes contributions from 56 Ni decay and subse- 
quent 56 Co decay. We use the analytic expression 
(jPinto fc Eastmanll2000l : iBersten et alj|2011l) 



,{t) = e m M m e- t/tNi +e co Mi 
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Fig. 9. — Bolometric luminosities and effective tempe ratu res cal- 
culated using the one-zone model described by eq. 1261 . Each 
curve is labeled with a M e j, M^;, and the total energy input. See 
text for further details. The general trend is that the energy injec- 
tion of the propeller mechanism creates a brighter, more quickly 
evolving supernova. 

where M^i is the mass of 56 Ni synthesized, = 3.9 x 
10 10 ergsg^ 1 s _ \ t m = 7.6 x 10 5 s, e Co = 6.8 x 
10 9 ergs g _1 s -1 , and t Co = 9.8 x 10 6 s. The factor 
of e _t /* T in equation (|26|) takes into account that the 
material eventually becomes optically thin to gamma- 
rays (although this factor only leads to small changes 
at late times). We have tested this simplified model 
against a wide range of nuclear-powered explosion calcu- 
lations (lEnsman fc Wooslevlll988t llwamoto et"aD 119981 : 
iDarbha et al.ll2010[ ). and found qualitatively good fits to 
the timescales and magnitudes of the peak luminosity. 

For Figure [5] we integrate equation (|26[) numerically, 
setting E int = E sn = 10 51 ergs at t = 0. Each curve 
is labeled by different values for M e j, Mm, and the to- 
tal energy input (initial supernova energy plus the pro- 
peller mechanism). The initial radius is 5Rq as ap- 
propriate for a compact Wolf-Rayet progenitor. We set 
TofF = (L/47rr 2 cr) 1 / 4 , where r = Vft and a is the Stefan- 
Boltzmann constant. This temperature is only accurate 
up to a time t « t T . The first two models (solid and long- 
dashed lines) explore the effect of a high input energy. 
The curve labeled with M cj = 10M Q and M Ni = O.7M 
dotted line) is repre sentative of a "hypernova" model 
Nomoto et alJ 120011 ). For this model we use an mi- 
tial radius of lOi?^, c onsistent with massive helium stars 
(Wooslcy et al. 1995). The model with merely 10 51 ergs 
(short-dashed line) is meant to be representative of a nor- 
mal Type Ib/c supernova. 

From these calculations, we find the general trend that 
additional energy injection from the propeller mecha- 
nism creates a brighter, more quickly evolving supernova 
(it will also cool faster and show optically thin features 
sooner). Within the framework we have described, it is 
not necessary that these events produce more 56 Ni than 



SUPERNOVA FALLBACK ONTO MAGNETARS 



9 



average. We therefore expect propeller-powered Type 
Ib/c supernovae to be associated with a wide range of 
peak luminosities, but to generically exhibit high veloci- 
ties of - (1 - 3) x 10 4 km s- 1 . 

5.3. High-Mass Envelopes 

If the envelope is more massive and has a hydrogen- 
rich composition, the lightcurve evolution can be sig- 
nificantly different, as is seen for Type IIP super- 
novae. The analytic featur es of th ese li ghtcurves 
were well summarized by iPopovl (|1993f ). whose 
work was confirmed an d expanded upon by the 
numerical simula t ions of lEastman et al.l (|1994l) and 
iKasen fc Wooslevl (|2009l) , and also studied with the first 
no n-LTE time-dependent r adiative-transfer simulations 
bv lDessart fc Hillierl (|2011[ ). The general picture is that 
the backward progression of a hydrogen recombination 
wave through the expanding ejecta causes the supernova 
to radiate at a fixed effective temperature set by the ion- 
ization temperature T c g = 2 1 / 4 Tj on . This continues until 
the entire envelope has become neutral, which truncates 
the lumin osity, r e vealin g 56 Co decay if it is sufficiently 
available. Popov (1993f) demonstrated that a hydrogen- 
rich envelope will exhibit a plateau phase when a cer- 
tain dimensionless parameter is greater than unity. We 
rewrite this condition in terms of a critical mass, finding 

M ei > 6El^ 5 R%^ A 3 { 3 T-% 3 M e , (29) 

where ko.34 = k/0.34 cm 2 g _1 is the electron-scattering 
opacity for a solar composition, Rq is the initial stellar 
radius with -R500 = i?o/500i?Q, and T5045 = Ti on /5045 K, 
corresponding to an effective temperature of 6000 K. For 
ejecta masses larger than this we expect a pro minent 
plate au phase and, scaling the analytic results of (jPopovl 
1993) to our values, with a plateau luminosity of 

L plat = 2.8 x 10 43 M- 1/2 ^ 6 5 ^/^S^T 5 4 / 4 3 5 ergs s"\ 

(30) 

and a plateau timescale 

Vat = ^M^E^rI^kI^T'I^ 3 days, (31) 

where we note that IKasen fc Wooslevl ()2009h find a 
slightly stronger scaling of t p \ a t oc _E -1 / 4 in their nu- 
merical results. The large energy input would also result 
in higher velocities of Vf ~ 10 4 km s -1 (scaling eq. [23] 
to a mass of ~ 10 — 2O-M0), which although not as high 
as in broad-lined supernovae, would be anomalously high 
for a Type IIP supernova. The high luminosities we find 
are similar to what is seen for many Type Iln super- 
novae (see Fig. 3 of Smith et al. 2008), but our events 
would not have nebular features from the interaction with 
winds and thus would appear distinct from Type Iln su- 
pernovae. 

To better demonstrate the impact of this energy injec- 
tion on the plateau phase, we plot examp le lightcurves 
using the analytic results of IPopovl (|1993| ) in Figure [T0l 
Beyond the plateau stage, the lightcurve may reveal a 
power-law decline from 56 Co decay, which we include for 
a range of 56 Ni masses, using equation ([2"5]). We do not 
plot the effective temperature since it is nearly constant 
at ~ 6000 K throughout the plateau phase. 



T 




Time (days) 

Fig. 10. — Bolomctric luminosity calculated using the analytic 
model of Popov (1993). The thick lines are models with an energy 
injection of 3 X 10 52 ergs, and the thin lines use 10 51 ergs (repre- 
sentative of a normal Type IIP supernova) . The additional energy 
injection results in a la rger l umi nosity, but shorter timescale, con- 
sistent with eqs. (130 P and (13 1 D - Beyond the plateau stage, the 
lightcurve is dominated b y 56 Co decay, which we plot for a range 
of 56 Ni masses using eq. (1 28 P as indicated (dotted lines). 

6. CONCLUSIONS AND DISCUSSION 

We presented a study of the effect of supernova fallback 
accretion onto newly born magnetars. The combination 
of spin, magnetic field, and fallback rate was used to cal- 
culate the time evolution of the magnetar spin, estimate 
how much material was accreted, and determine whether 
the magnetar can expel enough material via the propeller 
mechanism to prevent collapse to a black hole. Strong 
magnetic fields and short spin periods are generally more 
advantageous for hindering black hole formation (but as 
Figures [3] and 2] show, there are subtle changes to this 
picture depending on details of the time-dependent ac- 
cretion rate). Even in cases that avoid becoming black 
holes, ~ 0.3-Mq or more of supernova fallback material 
is accreted, so we expect magnetars formed in collaps- 
ing massive stars to be more massive than the canonical 
~ IAMq neutron star mass. As discussed in iJTJ there 
are at least three observed cases in our Galaxy of mag- 
netars associated with ~ 30 — 40M Q progenitors. The 
propeller mechanism suggests a natural connection for 
why neutron stars associated with massive progenitors 
should have magnetar-strength fields. 

Quickly spinning magnetars have been discussed as 
promising candidate systems for gravitational wave pro- 
duction via the time-changing quadrupole moment cre- 
ated by dynamical or secular instabilities (see Corsi & 
Meszaros 2009, Ott 2009, and references therein). We 
conclude that there are two main cases that may lead to 
the emission of gravitational waves when fallback accre- 
tion is important. In the first case, if the propeller mech- 
anism is active (typically £> > 5 x 10 14 G), the magnetar 
must begin with a sufficiently short spin period by the 
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process of cooling and contraction, as is found for some 
models explored by Ott et al. (2006, in particular, see 
their summary of /3- values in Table 4). It will then emit 
gravitational waves until it is spun down by accretion on 
a timescale of ~ 10 — 100 s. Since low accretion rates 
would extend the timescale for gravitational wave emis- 
sion, and these correspond to more energetic explosions, 
our model predicts that gravitational waves (if present) 
are most likely in the most energetic events that do not 
collapse to black holes. We note that in such cases a 
gamma-ray burst may be created directly by the magne- 
tar, as explored by Metzger et al. (2010, and references 
therein). In the second case, when accretion occurs di- 
rectly onto the magnetar surface, the magnetar is spun 
up sufficiently to emit gravitational waves. But as we 
discussed at the end of §5, this only proceeds until the 
magnetar collapses to a black hole after < 50 — 200 s 
(see eq. [18] ) . The formation of a black hole and its sub- 
sequent accretion may then power a gamma-ray burst, 
again predicting a possible correlation between gravi- 
tational wave emission and a powerful electromagnetic 
event. Note however that in this case the gravitational 
waves would precede any sort of launching of a relativistic 
jet (in contrast to the model of Piro & Pfahl 2007, which 
predicts gravitational waves coincident with the prompt 
gamma-ray emission, although both processes can occur 
in the same event). 

When a magnetar is in the propeller regime, the ex- 
pelled material collides with supernova e jecta, shock- 
heatin g it and energizing the supernova. iMaeda et al.1 
(2007) proposed that some ultraluminous supernovae 
may be explained by dipole emission from a rapidly spin- 
ning magnetar, which was worked out in detail by Kasen 
& Bildsten (2010) and Woosley (2010). We emphasize 
that our model is very different from theirs. In their case 
the magnetar directly powers the observed supernova lu- 



minosity. In our case the spin energy is injected earlier, 
creating a faster evolving supernova. We explored two 
regimes where this energy input may have a direct ob- 
servational consequence: (1) in the case of a low-mass 
(< 5Af ), hydrogen-deficient envelope, the additional 
energy gives rise to a broadlined Type Ib/c supernova, 
or potentially, a hypernova, depending on the amount 
of 56 Ni synthesized, and (2) in the case of a massive 
(> 10M Q ), hydrogen-rich envelope, we predict an event 
similar to a Type IIP supernova, although brighter and 
with higher velocities. 

Our predictions for the lightcurves of energetic super- 
novae are independent of the actual mechanism for inject- 
ing the energy, requiring only deposition at early times 
(<C td)- Therefore, independent of our specific model for 
how the energy is produced, Type IIP supernovae that 
have similar lightcurves as we demonstrate in §5.3, along 
with high velocities, indicate an exceptional amount of 
energy injection. Indeed, SN 2009kf has a luminosity 
that implies an explo sion energy of « 2 x 10 52 ergs 
(|Botticella et al.l [20101) . The databases of high-cadence 
transient surveys , such as the Palom ar Transient Factor y 
(|Law et all [20091 ) or Pan-STARRS ([Kaiser et al.ll2002f ). 
may reveal a larger population of Type IIP supernovae 
that, although not as extreme as SN 2009kj, may still 
require energy input beyond what is typically available 
for a supernova. 
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APPENDIX 

NEUTRINO-COOLED ACCRETION COLUMNS 

In cases where r c > r m > R, material is magnetically channeled before reaching the magnetar's surface. This 
is traditionally called an accretion column in the study of accreting, magnetized white dwarfs and neutron stars 
(|Frank et al.l[l99 2). For a dipole field, sin 2 9/r is constant, so that the path of the flow is described by the equation 
1/^m = sin 2 9/r (for simplicity we assume an aligned rotator). At a radius r from the magnetar, the material is 
squeezed into an area 

A{r) ks 7rr 2 sin 2 6 « Trr 2 (r/r m ). (Al) 
Assuming that the flow comes in at approximately free-fall, the velocity and density are 

/2GM\ 1/2 M 
^=(— ) • P0 = 2M n 

where the factor of two is because there are two poles. From this we estimate 

win = 1-8 x lO^M^rT/a 1 /2 cm s _ \ (A3) 
where r 12 = r/12 km. Combining equations ©, ([All) , and (jA2[) . the density is 

po = 1.6 x 10 8 ^Ml'Jr^^M^ g cm- 3 . (A4) 

The flow will go through a shock before reaching the stellar surface. This is checked by estimating the Mach number 
M of the flow, 
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where c s is the speed of sound, 7 is the adiabatic coefficient, and pb/p is the ratio of the magnetic pressure to the 
pressure of the gas (including ideal gas, radiation, and degeneracy contributions). For typical parameters, 

M = lA^MtTM% U r%* ^ , (A6) 

where we set 7 = 4/3 for a radiation dominated gas. In comparison, Zhang & Dai (2010) find that M. <C 1 because 
they assume much larger densities of ~ 10 12 g cm -3 , but from continuity this implies an infall velocity much less than 
freefall (v- ln ~ 10 5 cm s^ 1 ) contrary to our expectations. In addition, we estimate the mean free path for proton-proton 
collisions to be ~ 10~ 7 cm (using eq. [3.20] from Frank et al. 1992), much less than the width of the accretion column, 
so we expect the shock to be collisional. 

The flow can be broken into two main regions. The hrst is a supersonic flow starting from the edge of the magneto- 
sphere and then moving toward the magnetar pole. Then there is a shock, below which the subsonic flow settles onto 
the star. The jump conditions at the shock interface are 

Psh = 7p , p sh = ^p s hvfn, v sh = ^v in , (A7) 
for a strong shock with 7 = 4/3. Therefore the post-shock density is 

Psh = 1.1 x 10 9 ^ 7 M^ 7 r^ 2 M% 7 g cm- 3 , (A8) 

and 

1/4 



3p. 



'sh 



T s h : 

\ a j 

= 1.1 x lO^^X 3 / 28 ^ 8 ^ 28 K, (A9) 

is the post-shock temperature. 

The radiative diffusion timescale is approximately t ~ npr/c ~ 10 4 s, so the flow ca nnot cool via photon s. Instead 
we consider neutrino cooling via electron-positron pair annihilation, which is given by (jPopham et al.lll999t ) 

g P airs - 5 x 10 33 !* ergs cm" 3 s" 1 , (A10) 

where Tn — T/10 11 K. The timescale for this cooling is 

aT 4 

Cpairs ~~ ~ 

Qpairs 

= 5.2 x W-^ /7 MrT /2 ^u /8 Ml2 5/28 b. (All) 
It is also possible that Urea cooling is important, given by 

?Urca - 9 x 10 33 p 10 T 1 6 1 ergs cm" 3 s" 1 , (A12) 

for a composition of protons and neutrons (at these high temperature helium is photodisintegrated) . The timescale 
for this cooling is 

aT 4 



<?Urca 

= 6.3 x 10-V 3 6/7 Mi;f /14 4 7/4 MI 2 15/14 s. (A13) 
Urea cooling dominates when turca < £pair, implying an accretion rate 

M > 2.9 x 10~ 2 fi^M^rli 10 M s -1 . (A14) 

This is a rather high accretion rate in comparison to what we consider, so it is sufficient to focus on pair cooling. The 
height of the shock above the magnetar surface is 

H s h w impair = l.3^z /7 M^2 1/28 rll /8 MZ% 5/28 km. 

(A15) 

The shock therefore occurs at a radius of r s h = R + -ff s h- Note that the r on the righthand side of equation (|A15I) 
corresponds to r s h, so this equation is only accurate as long as H s h < r 8 h ~ R- 

When iJ s h > i?, we need to take into account adiabatic compression of material as it moves toward the magnetar 
pole. This gives a higher temperature at the magnetar surface in comparison to the shock radius, by an amount 
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Fig. 11. — Critical radii as a function of accretion rate for the problem of a neutrino-cooled accretion column. A given fluid element 
moves inward in radius (from top to bottom on the plot) at fixed accretion rate. In this way, one can read off what processes the fluid 
element experiences. When it reaches r m (dashed line), its motion is determined by the magnetic field. If r m > r c , then it will be expelled 
(we plot r c for 1 and 10 ms spin periods as examples; dotted lines). If r m < r c , then the flow will be channeled toward the magnetar pole 
and undergo a shock at r s h (solid line). It finally reaches the magnetar surface at R (dot- dashed line). As a comparison we also plot the 
approximation give by eq. J A 15 I t as the line labeled R + H^. 

T = T s h(i?/r s h) _1 , which revises the pair cooling timescale by a factor of (i?/r s h) 5 . We again write an equation for 
the shock height, 

H s u = r sh - R w W sh ipair(r s h)- (A16) 

This expression can be solved numerically for r s h, which we plot in Figure[TT]in comparison to other critical radii. This 
shows that for all accretion rates of interest, r s h < r m , so the location of the shock is always within where tunneled 
flow begins, as needed for consistency. 
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